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Abstract 

In this paper we review and refine a technique of Rioul to deter- 
mine the Holder regularity of a large class of symmetric subdivision 
schemes from the spectral radius of a single matrix. These schemes 
include those of Dubuc and Deslauriers, their dual versions, and more 
generally the pseudo-spline and dual pseudo-spline schemes. The tech- 
nique also applies to convex combinations of these, and therefore to 
^ several 'tension parameter' schemes. 
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^ 1 Introduction 

> 

Subdivision is a recursive method for generating curves, surfaces and other 



geometric objects. Rather than having a complete description of the object 
of interest at hand, subdivision generates the object by repeatedly refining 
its description starting from a coarse set of control points. Since subdivi- 
sion schemes are often easy to implement and very flexible, they provide a 
powerful tool for modelling geometry. However, analyzing their smoothness, 
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or regularity, can be difficult. The purpose of this paper is to review and 
refine a method proposed by Rioul [H] to determine the Holder regularity of 
a surprisingly large class of subdivision schemes from the spectral radius of 
a single matrix. A joint spectral radius analysis is not required. 
Consider the scheme 

fj+i,k — ^ &k-itf (1) 
t 

with finitely supported mask a = (a k ) k( zz, and coefficients a k G R, acting on 
the initial data fo± G M, k G Z. At each subdivision level j > 0, let fj be 
the piecewise linear function with value fj >k at the point 2~%. The scheme 
is then convergent if it has a pointwise limit / := lim^oo fj. We will assume 
that only a finite number of the initial data /o,fc are non-zero, in which case / 
has compact support. In the special case of the cardinal data / ,fc — $fc,o> the 
support of the limit / will be the interval [K, L] if clk, cll ^ and a k = for 
all k < K and k > L. We note that shifting the a k merely shifts /: if / is the 
limit of the scheme with shifted coefficients a k '■= dk+i then f(x) = f(x + l). 
The Laurent polynomial 

A{z) = Y,a k z k (2) 

k 

is the symbol of the scheme. We note that multiplying or dividing A(z) by 
any power of z results in the symbol of the scheme with shifted coefficients, 
and therefore has no effect on the regularity of /: if for example a k := a k+ \ 
then A[z) = z~ 1 A{z). It has been shown, for example by Dyn and Levin 
|13j . that a necessary condition for convergence of Q is that 

2j a 2k = 2j a ^k+l = 1, (3) 
k k 

and so we will make this assumption. This condition can be expressed in 
terms of the symbol as 

A(-l) = and A(l) = 2. (4) 

Let us now suppose, after shifting the coefficients a k as necessary, that 
A(z) admits the factorization 

A(z) = 2- r (l + z) r+1 B(z) (5) 
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for some r > 0, where b, the masking corresponding to B, is symmetric 
about bo, i.e., bk = fe-fe. Then the Fourier transform of b, 

b(f) := B{e~% e e R, (6) 

which is periodic with period 27r, is also real. 

Rioul showed in [18] that a /ou>er bound on the Holder regularity of (TO) 
can be determined from the spectral radius of a single matrix if b > 0, i.e., 
if b(£) > for all £ G [— 7r, 7r]. A surprisingly large class of schemes are 
of this type. For example, one can easily check that they include all the 
pseudo-spline and dual pseudo-spline schemes, from the explicit formulas for 
A(z) in H El El EJ. 

Rioul further showed that in the special case that the scheme is in- 
terpolatory, the lower bound is optimal. We will show more: that the lower 
bound is optimal whenever the cardinal function of the scheme has £°°-stable 
integer translates, which includes the interpolatory schemes as a special case. 
Using a characterization of such stability due to Jia and Micchelli [T7j, this 
leads us to conclude that the lower bound is optimal under the slightly stricter 
condition that b > 0. Such schemes include again all the pseudo-spline and 
dual pseudo-spline schemes. 

To be more specific, recall that a function such as / has Holder regularity 
a, < a < 1, written / G C a , if 

\f(x)-f(y)\<C\x-y\ a 

for all x,y el, and we write / G C q+a for q e N , < a < 1, if / G C q , i.e., 
/ is q times continuously differentiable, and g C a . Correspondingly, we 
shall say that the scheme ([T]) has Holder regularity 7 for some real 7 > 0, if, 
for < 7, / G C 13 for all initial data, and, for /3 > 7, / G" for some initial 
data. 

By the assumption of the symmetry of b, it must have the form 

b = (...,Q,b p ,...,bi,b o ,b 1 ,...,b p ,0,...), b p ^0, (7) 
for some p > 0, in which case 

p 

b(0=&o + 2^6 fc cos(A;£). (8) 
fc=i 
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Let Mb be the p x p matrix defined by 



if J = 0; 



b\i-2j\ + h+2j otherwise, 



(Mb)i+lJ+l - 

so that the first few examples of Mb, with p — 1, 2, 3, 4, are 



(9) 



N 



b 2b 2 



b 


2b 2 





h 


h + h 


63 


b 2 


bo 


b-2 



'bo 2b 2 26 4 0' 

h h+h h 

b 2 bo + 64 b 2 64 

63 bi bi b 3 



(10) 



Let p be the spectral radius of Mb- Our main theorem is 

Theorem 1. Ifh > and p > 1/2, a lower bound for the Holder regularity 
of ^ is r — log 2 (p). Ifh > this bound is optimal. 

In practice one can take the largest possible r in ([5]), so that the matrix 
Mb is relatively small and has relatively large spectral radius. As an example, 
consider the quintic Dubuc-Deslauriers scheme with mask 



256 



, 0, 3, 0, -25, 0, 150, 256, 150, 0, -25, 0, 3, 0, . . .)• (11) 



There is a factorization (|5]) up to r = 5, in which case one finds 



b = (6_ 2 , 6_i, 6 , h, 6 2 ) = -(3, -18, 38, -18, 3) 



and 



38 6 
-18 -It 



Writing s = sin 2 (£/2) yields 



b(f) = - (38 -36cos£ + 6cos2f) = 1 + 3s + 6s 2 > 



(12) 
(13) 

(14) 



for any ( 6 E. Since Mb has the two eigenvalues 9/2 and —2, Theorem [T] 
implies that this scheme has regularity 



5-log 2 (9/2) « 2.8301. 
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The remainder of the paper is organized as follows. In the next three 
sections, we recall some notation and results on the convergence and regu- 
larity of ([I]), and their relation to the growth rate of the 'derived scheme' 
defined by B(z). The following three sections explain how to estimate this 
growth rate in the case that b > and the section after that shows that this 
estimate is sharp if b > 0, proving Theorem [lj Finally, using the formulas 
of jH [9j [7] , we apply the theorem to compute the Holder regularity of the 
pseudo-spline and dual pseudo-spline subdivision schemes, and some convex 
combinations of these. 



2 Regularity 

The regularity of / is related to the behaviour of the divided differences 
of the scheme, and their differences, defined recursively as /j ] = f)^ and 

dfl = fj,k ~ fj,k-i, and for s > 1, 

fW — —J s ^ „W — fW _ fW 

Jj,k - g 9j,k ' 9j,k - Jj,k Jj,k-v \ 10 ) 

If A(z) is divisible by (1 + z) r+l for some r > 0, as in (JsJ) , then there is a 
scheme for the and g^ k for s = 0, 1, . . . , r. For such s, if we define their 
associated Laurent polynomials as 



k 



then 

and 

where 

and 



FlUz) = AW(z)Ff(z\ (16) 
df, 1 (z) = B»(z)df(z\ (17) 
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from which we obtain the derived schemes 

/l+i.fc = y^, a k-2£fgi ( 18 ) 

£ 

and 

= ( l9 ) 

Under this divisibility condition, it can be shown [T3] that a sufficient 
condition for / to be in C r is that 

\g%\ < c\ j , (20) 

for some constants C and A < 1, independent of j and k. Moreover, if 
1/2 < A < 1, then / e C r ~ lo ^ (x) . 



3 Reduction procedure 



How can we use (20) in the case that it holds with A > 1? Then we do not 
know whether / e C r , but if r > 1 we can use the 'reduction procedure' of 
Daubechies, Guskov, and Sweldens |2| to obtain information about lower or- 
der derivatives. Although the procedure was shown to work for interpolatory 
schemes in [2], it also applies to the more general scheme Q. 

Lemma 1. Suppose that A(z) is divisible by (1 + z) r+1 for some r > 1. If 



l[2(fy holds with A > 1 then 

\g [ ;, 1] \ < D 1 2^X\ (21) 



while if it holds with A = 1, 
for constants Di, D 2 , D 3 . 



g [ ;: 1] \<(D 2 + D 3j )2-\ (22) 



Proof. By the divisibility assumption, A^(— 1) = 0, and by the assumption 
that A(l) = 2 in Q, it also follows that A^(l) = 2. Therefore, 

E4! = E4! + i = i> (23) 
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and from (18) using summation by parts there is a constant C\ such that 
\flii,2 k+s -fll\<C^\g^\, 5 = 0,1. 

So, for any j > 1, if we represent any k G Z in binary form as A; = fcj, where 

kt = 2k t -i + s e , t = j, j - 1, . . . , 1, 
for some k G Z and Si, . . . , Sj G {0, 1}, then 

i/JS - 4ii < E \& - /i-i,^-J ^ ^ + A + • • • + A ^)- 

Hence, 

l/JSl <C 2 + C 1 C(1 + A + --- + A^ 1 ), 

and since 

this gives the result in the two cases A > 1 and A = 1. □ 



By applying this procedure recursively, it follows that if (20) holds for 
any A with 1/2 < A < 2 r then / G C r ~ log2X if log 2 A is not an integer, and 
/ G C r ~ log2 A_e for any small e > if log 2 A is an integer. 

4 Rioul's method 



k ■ 



With r in (5) now fixed, we have B(z) = B^'\z). If we also write bk = b 
9j,k — gfi, an d Gj(z) = the derived scheme for the can be expressed 



as 



9j+i,k = z2 b k-2e9j,e (24) 

t 

or in terms of Laurent polynomials as 

G J+l (z)=B(z)G J (z 2 ). (25) 

Rioul showed that if B(z) has the form ^ and if b > 0, (20) holds with A 
the spectral radius of a single matrix whose elements are formed from the 
coefficients bk- Note that if p = we must have bo — 1 and so we can take 
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A = 1. In fact in this case the scheme Q is the B-spline scheme of degree r 
and this merely confirms the well-known fact that the limit /, being a spline 
of degree r, belongs to C 13 for any (5 < r. Thus, we assume from now on that 
p > 1. Iterating (25) gives 

G J (z) = B J (z)G (z 2J ), (26) 

where 

Bj(z) := B(z)B(z 2 ) ■ ■ ■ Biz 2 '' 1 ). (27) 

But then 

B j+1 {z)=B{z)B j {z 2 ), (28) 
and so Bj(z) is the Laurent polynomial of the data bj^, where &o,fc = <5fc,o and 

h i+l,k = ^2 b k-2ibj,£. (29) 



In particular, bi^ = b^. Since ([26J) can be written as 

9j,k = bj,k-2iego,i, (30) 



it follows that 



\gj,k\ < max \b j>r \ 22 



9o,e\ 



and so (20) holds if there is some constant C such that 

max | < C'X 3 . (31) 

By induction on j, the values b^k are zero whenever k < —pj or k > pj, 
where pj := (2 J — l)p. 

Lemma 2 (Rioul). If B(z) in (5) has the form |?j) and b > then 
max | bj y k | = \bj t o\ for all j > 0. 

k ' ' 

Proof. Since 

Pj 
k=-pj 
8 



is a Fourier series, we have 



2tt 



Therefore, 



By the symmetry of the fefc, it follows from (29) by induction on j that the 
fej jfe are symmetric for all j, i.e., = bj t k- Therefore, bj is real, and by 
induction on j from (28), bj(£) > for all £. It follows that 



l fo 7,fc| < 



2tt 



b,(Ode 



□ 



5 A folded matrix 



Under the assumption of Lemma |2| it follows from (31) that (20) holds if 

\bj,o\ < C'X J (32) 

for some constant C Thus we only need to estimate the grow rate of the 
central coefficient bjo as j — > oo. Since bj_^ = bj^, one way to determine 
this rate is to study the behaviour of the vectors of coefficients 

bj = (b jfi ,b jtl , . . . ,b j:P _ 1 ) T , 



since they are self-generating and include bjQ. Indeed, from (29) 



bj+i,k — bkbj t o + ^^(&fc-2£ + bk+2i)bj,t, 

and using the fact that b-k = bk implies 

bj+i,k = bkbjfl + /~J(fc[fc-2i| + bk+2e)b 



and it follows that b J+1 = Mhj, where M is the p x p matrix Mb defined in 
©. Therefore, with e = (1, 0, . . . , 0) T G W, 



•j,o\ 



|M J e|L < \\M 3 \ 



p{M) 3 



as j — > oo. 



where p(M) is the spectral radius of M. Thus (20) holds with A = p(M). 
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6 An alternative folded matrix 



Rioul obtained the growth rate of bj$ in a different way, working from an 
alternative to equation (28). From (27) there is an alternative recursion: 



and therefore 



B j+l {z) = B{z 23 )B j {z) 



(33) 



This means that the subset of coefficients bj^ whose indices increase in steps 
of 2 J , rather than 1, i.e, Cj t k '■= bj^k, satisfy the equation 



c j+i,k — bgCj t 2k-e — bih-iCj,, 



(34) 



Since bj >0 = Cj >0 , we can obtain the growth rate of bj >0 by studying the be- 
haviour of the vectors 



instead of the hj. Similar to (29), equation (34) simplifies since the mask b 
is symmetric, for then Cj-i = Cj/, and so 

c j+i,k = b 2 kCjfl + '^^(b 2 k~e + b 2 k+e)cj/, 

and again because b-k = bk, 

c j+i,k = b 2 kCjfl + '^^(b\ 2 k-£\ + b 2k+ g)cj : e. 



Thus c 



j+i 



Ncj, where N is the p x p matrix defined by 



(N) 



>2t 



b\2i-j\ + b 2 i + j 



if J = 0; 
otherwise, 



(35) 



The first few examples of N, with p = 1, 2, 3, 4, are 
r, i \b 26i 



X 


26i 


2b 2 


b 2 


&i + h 


bo 





h 


b 2 



b 


26i 


26 2 


26 3 


b 2 


6l + &3 


6 + &4 


bi 


h 






bi 








6 4 


h 



10 



Therefore, 



\bj,o\ = \cj,o\ = ll^ e lloo = II^IU ~ p{NY, as j ->■ oo. 

Thus we can instead take A = p(N). However, because the transpose N T is 
similar to M, one has p(N) = p(M). 

7 Optimality 

Let denote the cardinal function of the scheme ([T|, i.e., its limit when the 
initial data is the cardinal data /o,fc — ^fc,o- Then the limit of the general 
scheme can be expressed as the linear combination 

Following Jia and Micchelli [17], we shall say that has £°° -stable integer 
translates if there is some constant K > such that for any sequence c = (ci)e 
in f°(Z), 

|| ^ c t <j)(- - ^)||l°°(r) > K\\c\\i°°(z)- (36) 

i 

Lemma 3. Suppose (ft has £°° -stable integer translates and f has regularity 
q + a for some q G No and < a < 1 . Then for any integer r > q, there is 
a constant C such that 

| 5 M| < C2^ r - q ~ a \ (37) 

Proof. As is well known, see e.g. the review by Dyn and Levin |13j . satisfies 
the two-scale difference equation 

<P(x) = J2 a k<P{2x-k), (38) 

k 

and therefore, for any j > 0, 

/(,■) ^/, / ,,(2',--/). (39) 

e 

We can use this equation to relate any divided difference of / of the form 
/g := [2-'(y - q), 2~\y - q + 1), . . . , 2^(y)]f, 
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for y G R, to the divided differences of the scheme. Putting x = 2 J — k) 
in (39) gives 

f{2- i (y-k))=Y,fj+-h<Kv-Q, 

i 

and, using the cases k — 0, 1, . . . , q, and the linearity of divided differences, 



Similarly, if 
then 



pcM ._ _ f[q] 



Recalling that / has compact support, if / has regularity q + a, there is 
some C > such that for any £ 0) £i e R> 

i/ (g) (6)-/ (9) (eo)i<ci6-eor, 

and, by a standard property of divided differences, for each j and y, 

i5Si = i/ (9) (6)-/ (9) (eo)iM 

for £o;£i ^ (2~ J (y — q — 1), 2~ J (?/)). Therefore, for any y, 



where 
Therefore, 



C" = C(g + l)7g!. 



and by (36) it follows that for any £ 6 Z, 



Finally, by applying the divided difference definitions (15) recursively, we 
obtain (37). □ 
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It follows from this lemma that as long as has £°°-stable integer trans- 
lates, the lower bound on the Holder regularity of the scheme that we ob- 
tained from the spectral radius p = p(M) in Section [5] is optimal. To see 
this, suppose that p corresponds to a real eigenvalue of M. After taking the 
initial data 

r] [r] \T 



(a 11 a 1 



in the derived scheme (|24|) as the corresponding eigenvector, extending by 

9olk = 9oi y iel( 

the rate p* , i.e., 



9o,-k = 9o \ yields a sequence of divided differences g^ k that grows at precisely 



\9V k \~P>. (40) 
In the other case that p corresponds to a pair of complex conjugate eigenval- 



ues, we again obtain the growth rate (40) by taking the real part of one of 



the corresponding eigenvectors. Suppose now that / G C r lo ^p+ e f or some 
small e > 0. Writing the exponent as 

r - log 2 p + e = q + a, 

for q G No and < a < 1, if p > 1/2, we have q < r, and the lemma then 
implies that 

I flfil < C2 i(log2 "- e) = Cp j 2- je , 



which contradicts (40). 



It remains to consider when has £°°-stable integer translates. If the 
scheme ([T]) is interpolatory, in the sense that a,2k = <5fc,0; then <f)(k) = 5k$ 



and so the stability condition (36) holds with K = 1. Thus, for the Dubuc- 
Deslauriers schemes the spectral radius of M yields the optimal regularity. 
Rioul showed this and tabulated the regularities of these schemes in [T8] . 

To show stability in the general case, we apply some results by Jia and 
Micchelli [T7j. We denote the (continuous) Fourier transform of <j) by 



0(0= / tWe-^dx, £e 
Jr 

Since the scheme ([!]) has constant precision, 

^2<p(x-£) = 1, xeR, 
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and, as shown by Jia and Micchelli [TTl Theorem 2.4], 0(0) = 1. Since the 
Fourier transform of (38) is 



0(0 = 2- 1 a(£/2)0(£/2), 

it follows that 

oo 

?(e) = n( 2_1 ^/ 2i ))- 

3=1 

A sufficient condition [TTJ Theorem 3.5] for to have £°°-stable integer trans- 
lates is that 

sup |0(£ + 2vrf) | > 0, for all (6i (41) 



Consider then again the case that the derived scheme (24) holds. Then 

a(0 = 2cos r+1 (e/2)b(0, 

where, since a(0) = 2 under the assumption of convergence, b(0) = 1. In the 
B-spline scheme of degree r we have b(z) = 1, in which case we can write the 
symbol as A r (z) = (1 + z) r+1 /2 r . The cardinal function r is the B-spline of 
degree r centred at 0, and we have 



. K)= ncos-'(^«)=(^) r+1 



It then follows that 

oo 

?(fl = ?r(flIIb(e/2'). 

3=1 



Since the condition (41) holds for the B-spline r we deduce that has 
stable integer translates if b(£) > for all (el. This completes the proof 
of Theorem Q] 

This idea of comparing with B-splines, which are known to be stable, was 
used by Dong and Shen [H Lemma 2.2] to show that all pseudo-splines are 
stable. 
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8 Pseudo-splines and dual pseudo-splines 



It turns out that a surprisingly large class of subdivision schemes that appear 
in the literature satisfy the conditions of Theorem [TJ including the pseudo- 
spline and dual pseudo-spline schemes. 

Consider the family of primal schemes with symbols 




It is shown by Dong and Shen [8] that these are the symbols of the pseudo- 
splines of type II, which were introduced by Daubechies, Han, Ron, and Shen 
[I] and studied systematically by Dong and Shen j9]. The 2m-point Dubuc- 
Deslauriers schemes form the subfamily with symbols A m)TO _i. Other choices 
of I and m lead to other known schemes. Dividing A m j by 2a m gives us the 
symbol for the shifted derived scheme, 

B m M = = £ ( m " 1 + k )5\z) = £ W, (43) 

k=0 ^ ' i=-l 

with 

A family of schemes analogous to the Dubuc-Deslauriers schemes was 
introduced by Dyn, Floater, and Hormann [11] and generalized by Dyn, 
Hormann, Sabin, and Shen [7] to a family of dual schemes with symbols 

A^(z) = ^a m (z)^( m ~ 1 ^ + k )8 k (z), m>l, I > 0. (45) 

fc=0 ^ ' 

Dividing A m ^z) by - J ^a m (z) gives us the symbol for the shifted derived 
scheme, 

BmM = B [ %\z) = j:( m ~ 1/2 + *W) = J2 b mM z\ (46) 

fc=0 ^ ' i=-l 
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with 

i 

b m ,i,i = {-iy 

For any m > 1 and I > 0, since 



m - 1/2 + fc 



2fc 
i + k 



,-2& 



<f(e*) = sin 2 (£/2) > 0, £ G [-tt, 

a. 



71 



we see that h m) i > and b mj ; > and so Theorem |l| can be applied. 

Example 1. Consider the primal schemes with 1 = 1. From (J^44) we 
find that M is a 1 x 1 matrix with spectral radius & m ,i,o = 1 + m/2 > 1/2. 
We conclude that the primal scheme with symbol A m> i has Holder regularity 
2m — log 2 (?^ + 2), which was also shown by Hormann and Sabin FR 



Example 2. For the primal scheme with (m, I) = (4, 3), one finds B^ 3 (z) 



+ b^z 3 , with 
(b , /;:>,) 



16 



;-5, 40, -131, 208, -131, 40, -5). 



The regularity of the original scheme with symbol is then 7 — log 2 (p) ~ 
3.55113, with p w 10.91976 the largest zero (in absolute value) of the polyno- 
mial 



det(M-AJ) = det 



b - A 
bi h 
b 2 



2b 2 
bo 





h 
bo- A 



217 

-A 3 + 7A 2 + A -125. 
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Similarly, one can, for any integers m > 1 and I > 0, compute the Holder 
regularities of the schemes with symbols A m j and A m j_ as the log 2 of algebraic 
numbers of degree at most I. For any integers 1 < / < m < 8, Tables [T] and [2] 
contain Holder regularities for the primal and dual schemes, showing that 
the first four decimals of the lower bounds presented by Dong, Dyn, and 
Hormann [6] are, in fact, significant. 

Moreover, since 



a(e^) = cos 2 (£/2) > 0, £ G 



-7T, 7T 



any convex combination of the schemes v4 m ;, i.e., the schemes with symbols 
(1 - w)A mi)h (z) + wA m2h (z), w e [0, 1], 
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m = 2 
m = 3 
m = 4 
m = 5 
m = 6 
m = 7 



I = 1 
2~ 

3.67807 
5.41504 
7.19265 
9 

10.83007 



I = 2 

2.83007 
4.34379 
5.92502 
7.55781 
9.23111 



I = 6 



I = 7 



3.55113 
4.96207 
6.43997 
7.97187 



4.19357 
5.53250 
6.93577 



4.77675 
6.06273 



5.31732 



m = 8 12.67807 10.93702 9.54804 8.39272 7.41006 6.56398 5.82944 



Table 1: Holder regularities rounded to five decimals for the primal schemes 
A mi i with 1 < I < m < 8. 







/ = 1 


I = 2 


I = 3 


I = 4 


I = 5 


I = 6 


I = 7 


m 


= 2 


2.83007 














m 


= 3 


4.54057 


3.57723 












m 


= 4 


6.29956 


5.12711 


4.24726 










m 


= 5 


8.09311 


6.73575 


5.69355 


4.85423 








m 


= 6 


9.91254 


8.38994 


7.19984 


6.22682 


5.41143 






m 


= 7 


11.75207 


10.08039 


8.75493 


7.65811 


6.72934 


5.93283 




rn 


= 8 


13.60768 


11.80033 


10.35034 


9.13861 


8.10385 


7.20968 


6.43070 



Table 2: Holder regularities rounded to five decimals for the dual schemes 
A m)[ with 1 < I < m < 8. 



have a derived scheme satisfying the conditions of Theorem [T} The same is 
true for the dual schemes A m j. 

Figure [TJ for instance, shows Holder regularities of convex combinations 
of the primal and dual Dubuc-Deslauriers schemes. Similarly, as remarked 
by Dyn et al. [7], for L > 1, the symbols of the schemes S 2 l(u!) of Choi, 
Lee, Lee, and Yoon pQ are affine combinations of Al,l-x and Al+i,l-i, while 
the symbols of their schemes S^L-xiw) are affine combinations of the duals 
Al,l-i and Al-i,l-i- When the combinations are convex, Theorem [l] can be 
applied. 
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